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NONLINEAR EQUATIONS FOR BENDING OF ROTATING BEAMS WITH 


APPLICATION TO LINEAR FLAP-LAG STABILITY OF 
HINGELESS ROTORS 


Dewey H. Hodges and Robert A. Ormiston 

Ames Research Center 
and 

U. S. Army Air Mobility R & D Laboratory 
Moffett Field, California 94035 


SUMMARY 


The nonlinear partial differential equations for the flapping and lead-lag degrees of freedom of 
a torsionally rigid, rotating cantilevered beam are derived. These equations are linearized about an 
equilibrium condition to study the flap-lag stability characteristics of hingeless helicopter rotor 
blades with zero twist and uniform mass and stiffness in the hovering flight condition. The results 
indicate that these configurations are stable because the effect of elastic coupling more than com- 
pensates for the destabilizing flap-lag Coriolis and aerodynamic coupling. The effect of higher bend- 
ing modes on the lead-lag damping was found to be small and the common, centrally hinged, spring 
restrained, rigid blade approximation for elastic rotor blades was shown to be reasonably satisfactory 
for determining flap-lag stability. The effect of pre-cone was generally stabilizing and the effects of 
rotary inertia were negligible. 


INTRODUCTION 


Helicopter rotor blades that are cantilevered to the rotor hub can substantially improve 
helicopter flying qualities by increasing rotor control power and angular damping. This type of rotor 
is known as a hingeless rotor because it lacks the flapping and lead-lag hinges found on conventional 
articulated rotors. The absence of such hinges requires the blades to bend elastically in two direc- 
tions: perpendicular (flapping) and parallel (lead-lag) to the plane of the rotor for centrifugal relief 
of aerodynamic forces. In general, the bending deflections are also accompanied by torsional 
deformations of the blade. 

Although experience with hingeless rotors is not extensive, previous research has indicated that 
for certain configurations, the coupling between the flapping and lead-lag degrees of freedom may 
result in unstable oscillations (refs. 1 and 2). This condition is known as flap-lag instability. It is a 
restricted case of the general hingeless rotor blade stability problem because the torsion degree of 
freedom is not included. This restricted problem is quite useful, however, because it gives consider- 
able insight into the coupling mechanisms of flap-lag motion which can contribute to rotor blade 
instabilities. These coupling mechanisms are primarily due to the combination of aerodynamic 



and inertial coupling (Coriolis and centrifugal), which is essentially proportional to the mean flapping 
or coning of the rotor blade. Since the coning increases with rotor thrust, flap-lag stability is of most 
interest for moderate to high thrust conditions. The Coriolis and centrifugal inertial forces are basi- 
cally nonlinear but can be linearized for small oscillations about the mean flapping and lead-lag 
deflections. 

The prominent feature of previous research (refs. 1 and 2) is a method of approximating the 
continuous elastic bending deflections of a hingeless rotor blade with a centrally hinged, spring 
restrained, rigid blade. Such an approximation reduces the problem from one of partial differential 
equations to one involving only ordinary differential equations. The solutions then become simpler 
to obtain and the results are somewhat easier to interpret in physical terms. Another important 
reason for using the approximate rigid blade representation is that equations suitable for the elastic 
bending of rotating beams including the nonlinear effects noted above have not been available, 
whereas the nonlinear rigid blade equations are relatively easy to derive. The most general linear 
treatment for the elastic rotating beam may be found in reference 3. 

The purpose of the present report, therefore, is to derive a suitable system of nonlinear 
equations for the flapping and lead-lag bending degrees of freedom of a torsionally rigid, rotating 
cantilever beam. A general treatment of the nonlinear inertial forces requires an additional degree 
of freedom for radial deflections of the blade and use of a nonlinear radial strain equation valid for 
large bending deflections. The radial degree of freedom is ultimately expressed in terms of the two 
bending deflections so that the final equations contain only the two primary degrees of freedom. A 
procedure for linearizing the equations about a particular equilibrium condition is given and results 
for the stability of various untwisted rotor blade configurations with uniform mass and stiffness are 
presented in the hovering flight condition. In addition to revealing the stability characteristics of 
elastic hingeless rotor blades, these results will also serve to establish the accuracy of the common 
rigid blade representation of a hingeless rotor blade. 

One result of reference 2 indicated that flap-lag instabilities were highly sensitive to the elastic 
coupling between the flap and lead-lag degrees of freedom. For the rigid blade representation, this 
coupling could be controlled by varying the orientation of the flap and lead-lag restraint springs at 
the blade hinges. For the actual elastic blade, the analogous elastic coupling effect is dependent on 
the orientation of the elastic principal axes of the rotor blade cross section and the radial distribution 
of bending stiffness. The results of reference 2 showed that when the restraint springs of the rigid 
blade model were configured to represent untwisted, uniform stiffness cantilever blades, the flap-lag 
instability was suppressed. Since the present numerical results are restricted to uniform untwisted 
blades, it is expected that unstable oscillations will not be encountered. However, in addition to 
flap-lag stability, the damping of the lead-lag bending motion is of considerable practical importance 
because the inherent aerodynamic and structural lead-lag damping of hingeless rotors is very small. 
This leads to potentially high fatigue stresses and vibration, and increases the susceptibility of the 
helicopter to coupled rotor-airframe instabilities. The results to be presented will compare lead-lag 
damping of the rigid blade representation with the elastic blade and examine the sensitivity of the 
elastic blade damping to the number of bending modes included. 
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EQUATIONS OF MOTION 


The equations of motion for a rotating cantilevered beam, rigid in torsion, are derived using 
Hamilton’s variational principle. They include three nonlinear equations for radial, lead-lag, and flap- 
ping displacements, u, v, and w, respectively. Aerodynamic forces are represented as a distributed 
generalized loading and are derived from simple quasi-steady strip theory. The initial derivation is 
valid for beams with variable mass and stiffness properties, pre-cone, and twist; however, for the pres- 
ent purposes it is sufficient to consider only untwisted configurations with uniform mass and stiff- 
ness. The resulting nonlinear equations can be linearized for small perturbations about a steady-state 
operating condition. This enables the radial displacement u to be eliminated from the equations of 
motion leaving only the primary flap and lead-lag displacements. The final perturbation equations 
are linear variable coefficient integro-partial differential equations, which are solved by Galerkin’s 
method. 


Variational Formulation 


Coordinate systems— The X, Y, Z coordinate system, figure 1, rotates with constant angular 
velocity about the Z axis and its origin is fixed in inertial space. The plane of rotation is defined by 
the X, Y axes. The position of a beam S 
is defined in a second x, y, z coordinate 
system, with the x axis coincident with 
the undeformed position of the beam and 
rotated about the Y axis through an angle 
/3 pc , the pre-cone angle. 

The displacements of the beam 
elastic axis u, v, w are defined along the 
unit vectors n i , n 2 , n 3 , which are parallel 
tox, y, z, respectively. The beam is canti- 
levered at the origin and the elastic 
axis, aerodynamic centers, and mass cen- 
ters are coincident with the x axis in the 
undeformed state. The beam pitch angle 
9, which defines the orientation of the 
cross section principal axes y' and z', is 
shown in figure 2. For an arbitrary dis- 
placement, the position vector of a point 
on the beam will be 


> / dv 9w\_> 

= [ x+u ~ y ^~ z xh 


+ 0+y)n 2 + (tv + z)n 3 



( 1 ) 



Figure 2 — Blade cross-section coordinate system. 
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This equation assumes that squares of bending slopes are negligible with respect to unity and 
that plane sections remain plane. The radial elongation of the beam 9w/9x is assumed to be of order 
(bv/bx) 2 and ( bw/bx ) 2 and thus can also be neglected with respect to unity. Similar reasoning applies 
to terms of order (y/R ) 2 and ( z/R ) 2 . 


Kinetic energy— The velocity of a point on the blade with respect to inertial space is 


.♦6 ?•->■_> 
v = — + X r 
bt 


( 2 ) 


The first term is the velocity with respect to the rotating x, y, z coordinate system. 


b7 _ / bu 
bt~ \bt~ 


b 2 v b 2 w 

y Z 

9x bt 9x bt 


bv bw 
In, + — n 2 + — n, 
) 1 bt 2 bt ' 


The second term due to the coordinate system rotation is 


nxr = 


f / bv bw\ 

= -f2cos P pc (v+y)n 1 + cos j3 pc fx + u -y — -z — I 

- n sin (3 pc (w + z ) ] n 2 + sin j3 pc (v + y)n 3 


The velocity is therefore 


[ bu 

a7' 


b 2 v b 2 w 

n cos &p C (v + y) -y r — — - z 


dx bt bx bt 




[ bv ( bv ___9w\ 1 ^ 

— + n cos /3 p Jx + u - y — - z — J-n sin 0 pc (w+z) l n 2 


3w , , 

+ — +f2sin(3 pc (v + y) Jn 3 




The kinetic energy is 


-m 


pv • v dx dy dz 


(3) 


(4) 


(5) 


( 6 ) 


4 



and its first variation is 


8T 


- Iff 


pv • 8v dx dy dz 


(7) 


where 


-► (88u d 2 8v d 2 8w 

5v = I f2 cos (3 nr 8v -y - — - — z - — — 

\ 9 1 P c dx dt dx dt 




f 35v / 35v 35w\ *]_> 

L a7 +ncos V 5w J n2 


( 


asw . 

+ I + S2 sin fi nr Sv 

a t p c 


)«, 


( 8 ) 


From the cross section geometry of figure 2 

y - y cos 6 -z' sin 9 
z = y sin 6 + z' cos 9 


(9) 


Substituting equations (5), (8), and (9) into (7) yields the following expression for 8T. All 
quantities of order (dv/dx) 2 are ignored with respect to one. 

/ R i/du \ d8u (dv _ \35v 

m | \di' nv cos ^p c ) 17 + la7 +nx cos Ppc + nucos Ppc' nwsin Ppc Jjy 

0 

( a w \ dsw r 3v "j 

— + £2v sin P pc J -jj- + — + i2(x + «)cos p pc - fiw sin (i pc cos p pc 8u 

+ [(i7 + nv sin sin ~ (l7 ~ cos ^p c ) n cos v] 5v 


^ + fl(x + m)cos P pc - £2w sin j3p C J tt sin P pc 8w | dx 


( 10 ) 


where 


m 


= J jp dy' dz' 


(ID 
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Note that due to symmetry properties of the cross section, the following integrals are equal to 

zero 

ffy' d y ' d z\ ff z ' dy' dz', fj’ 'z' dy' dz' 

Because higher order terms were discarded, the kinetic energy expression does not include 
contributions of rotary inertia. Although this is a very small effect for rotor blade configurations, the 
rotary inertia terms were included for calculation of certain results to be discussed below. 

Strain energy — Previous derivations of rotor blade equations restricted to the case of small 
deflections have employed linear strain-displacement relations. Because the present derivation is 
addressed to the flap-lag stability problem with moderately large displacements, the linear relations 
are not valid. The following strain-displacement relation is used, which has been modified from 
reference 3 by the inclusion of the underlined nonlinear terms: 


3 u 1 (bv\ 1 /dw\ ,/3 2 v 3 2 w 

= — +-(— 1 +-( — ) -y — - cos 6 + — - 

3x 2 \3x / 2 \ 3 xj \bx 2 bx 2 



( 12 ) 


For long, slender beams, the strain energy is independent of the shear modulus and is simply 


U 


= fff^r ixdy,,z 

S 


(13) 


The first variation is 


8U = j J ~ J" Eebe dx dy dz 
5 


(14) 


where 


35w 3v 35v 

8e = -i + 

bx bx bx 


3w 35w ,/ 

bx bx ^ \ 


3 2 5v 3 2 8w 

_ cos9+ _ stafl 


) +z (' 


3 2 5v 3 2 8w \ 

_.^ S ._ cos9 j (15 ) 


Substituting equations (12) and (15) into ( 1 4) and using equation (9), we obtain : 
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R 


8U = 


-J 


' EA 


3 u 1 / 9v \ 2 1 / 3w\ 

dx 2 1 3 x J 2 1 3 x J 


'd8u + 3v 3Sv 3w 36w 
V 3x 3x dx dx dx , 


( d 2 v 3 2 w \(d 2 8v 3 2 5w ' 

+ El 2 ( r cos 6 + — r- sin 0 — — cos d + ■ ■ ■, sin 6 

\ dx 2 dx 2 / \ 3x 2 3x 2 j 


( d 2 v 

t£,1 (a? si 


3 2 w \(d 2 8v d 2 8w 

sin 6 - — - cos d - - sin d — — 7 - cos 6 \(dx 
dx 2 / \ 3x 2 dx 2 


( 16 ) 


Nonconservative forces— The virtual work 8W of the nonconservative forces is usually expressed 
as a time integral: 


t 2 1 2 

J' 8W dt = J" Q • 8q dt 


(17) 


where (5 is the vector of applied distributed loads and 8~q is the vector of virtual displacements. 
Hence, 


to R 


f -If {L x 8u + Ly8v + L z 8w)dx dt 


(18) 


\ 0 


where L x , Ly, L z , and 8u, 8v, 8w are the components of Q and 8q, respectively. 


General equations— Hamilton’s variational principle is stated as follows: 

[8T- 8U+ 8W]dt = 0 


/ 


(19) 


We now substitute equations (10), (16), and (18) into equation (19) and perform, by parts, all 
integrations that are necessary to obtain an expression of the form 



h 0 


)8u + ( )5v + ( )5w] dx dt + boundary terms = 0 


( 20 ) 
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For arbitrary, admissible variations 5 u, 5v, 5w, the three expressions in parentheses must vanish 
individually (as well as the boundary conditions). This entire series of operations yields the following 
three nonlinear partial differential equations: 


9 / 9v „ \ 

" Sx iEAe) ~ m 2n a7 + n 


(21a) 


9 

/ 9v ^ 

9 

, , 9 2 v 

9 2 vv 

9x 

[EAe — 

\ dx / 

+ 9x 2 

(El 2 cos 2 0 + El . sin 2 0 ) — - 

dx 2 

+ (El 2 - EI t )sin 0 cos 0 — — 

dx 2 


+ m 


9 2 v 

— r- - fi 2 v + 2fi 
dt 2 


/9u 

\9? _l 


' pc 



(21b) 


9 

/ 9w\ 

9 

9 2 v 

\ 

£ 

9x 

( ) 

9x 2 

(El 2 - E^ )cos 0 sin 0 — - 

9x 2 

+ (EI 2 sin 2 0 + Eli cos 2 0) — — 

9x 2 


(d 2 w dv\ . 

* m \w + n V 1 + 2n, V a? ) * <21c) 


where 


9m 1 / 9v \ l/9 w 

9x 2 \9x y 2 y9x 


(2 Id) 


Note that e is the strain at the elastic axis. In obtaining the equations, the small angle assumption 
has been invoked for the pre-cone angle, sin (3p C = (3 p C , @ 2 C « 1. These equations are nonlinear, 
partial differential equations of fourth order with variable coefficients in x, and of second order in t. 


Linearization of the Equations 

Before the equations are linearized, they are restricted to configurations with untwisted blades 
and uniform mass and stiffness properties. The equations will also be made dimensionless by dividing 
through by mfi 2 /?. The dimensionless displacements, u, v, w are then based on the rotor blade radius 
R, and the independent variables become £ = x/R and r = fit. The radial aerodynamic force compo- 
nent L x is of second order for small 0 and w', and consistent with second-order structural terms, can 
be discarded. The equations are: 
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(22a) 


de dv 

-K £ - 2 — =0 

3£ 3 t 



3 4 v 

ai 4 


+ Aj 2 


3 4 w 

ai 4 ” 


3 2 v _ /3 m 9w\ - 

dr 2 \dr P c 3 t y 


(22b) 



and 


3 u 1 /dv \ 1 / 3w\ 

ai + 2 \ ai ) + 2yaF / 

A, 2 = Aj 2 cos 2 0 + A 2 2 sin 2 9 
A 2 2 = Aj 2 sin 2 0 + A 2 2 cos 2 0 
Ai 2 = (A 2 2 - A, 2 )cos 0 sin 0 
A t 2 sEIt/m&R* 

A 2 2 =EI 2 /mn 2 R 4 
/sf = £4/mfi 2 2? 2 

/>j/ = Ly/fflEl 2 R 

Lz = Lz/mSl 2 R 


(22d) 


(22e) 


The terms defined in equation (22e) are the dimensionless aerodynamic forces and elastic 
section properties. Note that 0 is now the blade pitch angle, and is constant over the radius. The 
parameter K is the ratio of longitudinal elastic forces to centrifugal forces, whereas Aj 2 and A 2 2 
are the ratios of the section bending stiffness to the centrifugal forces. Examination of the equa- 
tions reveals that tensile forces are associated with K, while A] 2 , A 2 2 identify the elastic bending 
terms. These equations differ from previous results (ref. 3) in two ways. First, the radial force equa- 
tion includes the centrifugal term due to lead-lag velocity (-2(3v/aT)). This means that the tension 
(or strain e) terms in the flap and lead-lag equations are nonlinear, because e is a function of lead-lag 
displacement. Second, the retention of the radial displacement u in the lead-lag equation produces 
the Coriolis forces that are responsible for the destabilizing effects on flap-lag oscillations. The radial 
displacement u is related to the longitudinal strain e and the flap and lead-lag displacements by the 
nonlinear strain-displacement relation. 
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The nonlinear equations may be linearized about a suitable equilibrium operating condition 
without losing the essential features of the nonlinear flap-lag inertial coupling. The displacements are 
assumed to consist of the sum of time-independent parts and the time-dependent parts. These are, 
respectively, the steady-state equilibrium displacements and the perturbation displacements: 


t) = u 0 (Z) + Am(£, r) 
»>(£, t)=v 0 (£) + Av(£,t) 
w(£, r) = vv 0 (£) + Aw(£, r) 


(23) 


Also, it is assumed that 


e(£, t) = e 0 (|) + Ae(£, r) 
L/Z, T) = L yo tt) + AL y U, r) 
L z ( t, r) = r zo (£) + Al z (£, r) 


\ 


(24) 


/ 


The equations that result after substituting equations (23) and (24) into equations (22), and 
setting all perturbation quantities equal to zero, are called equilibrium equations. The equations that 
result from making the full substitution of equations (23) and (24), subtracting the equilibrium equa- 
tions and discarding higher order products of perturbation quantities, are called the perturbation 
equations. The equilibrium equations are 
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and the perturbation equations are 


3Ae 3Av 

-K 2 — = 0 

3£ 3t 


(26a) 


3 

-K — 


( dv 0 \ 3 / 3Av\ 

( Ae w )~ k k \° IF/ 


— „ 3 4 Av — 3 4 Aw 3 2 Av 
2 


+ A 2 2 t TF + A i'> — — + — — " Av +2 


3^ 


3£ 4 3r 2 


3 Am 

3 T 


3 

-a: - 


( dw 0 \ 3 / 3Aw\ _ 

e “dF/' A ai\ 6 ° “aF ) + Al 


3Aw _ 

-2P pc —=AL y (26b) 


— 3 4 Av - 3 4 Aw 3 2 Aw 3Av — 

2 — 5 - +A i r- + — =- + 2j3 n „ = AX_ (26c) 

2 3£ 4 1 3£ 4 3 r 2 P c 3r z 


where 


Ae = 


3Am dv 0 

aF + d? 


3Av dw 0 3Aw 

3T + dF ”aT 


(26d) 


Except for the equilibrium nonlinear strain equation (25d), the radial displacement u 0 does 
not appear in either set of equations. Therefore, a solution for u 0 is not required and equation (25d) 
is superfluous. The radial force equilibrium equation may be integrated and used to eliminate e 0 in 
the remaining equations. Equation (25a) yields 


e 0 


l-l 2 
2 K 


(27) 


Substituting this equation into equations (25b) and (25c) yields the following equations for flap 
and lead-lag equilibrium: 
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The perturbation equations may also be reduced to two equations by making substitutions for 
the terms e 0 , Ae, A m, and 9 A m/9t. The perturbation strain is evaluated by integrating the perturba- 
tion radial force equation (26a). 



e 


(29) 


Combining this equation with equation (26d) gives an expression for the perturbation radial 
displacement A u in terms of flap and lead-lag displacements 


9 Am 

~3F 


dv 0 9Av 

dF “9f 


dw 0 9 Aw 2 r l 9Av 

~d? 1 T + KJ 


(30) 


The final term may be omitted because K is typically on the order of 1 0 2 . This approximation 
is equivalent to the assumption that the beam is an inextensible member in terms of perturbation 
deflections. Integrating equation (30) then yields the perturbation radial displacement: 


A _ f^(dv 0 9Av dw 0 9Aw\ 

Am — / [ — H Id? 

J \d£ 9£ d£ 9£ / 


(31) 


and 


9Am 

9r 


9 r^/dv 0 9Av dw 0 9 Aw \ 

~VtJ y"dF "9F +_ dF ~df~ J ' 


(32) 


Equations (27), (29), and (32) may be substituted into equations (26b) and (26c) to obtain the 
final two perturbation equations in terms of the flap and lead-lag displacements Aw and Av. 


dAv /l-£ 2 \ 9 2 Av dv 0 9Av d 2 v 0 f 9Av — , 9 4 Av _ 9 4 Aw 9 2 Av 

IF T~ ) ~aF ^ 2 1? ~ 2 WJ — d£ + Az ~ + A ” — 


9£ 4 


2 9£ 4 9r 2 


- Av - 2/3 


9Aw 


P c 9r 


9 1 f 

9 tJ ^d£ 


9Av dw 0 9Aw\ 

+ — ld£ = A L v 

9£ d| 9£ » y 


)■ 


(33a) 
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1 


3Aw 




3 2 Aw dw 0 3Av 

r +2 — - 

3£ 2 d| 3 r 



3Av _ 3 4 Av T , 

ir d{+s , s i?- tS ‘ 


3 4 Aw 

3£ 4 


+ 


3 2 Av? 

3 T 2 


+ W pc 


3Av 
3 r 



(33b) 


The singly underlined portions of equations (33) represent the linearized Coriolis and centrifugal 
terms that result from retaining the nonlinearities noted above. They are necessary for a proper 
description of rotor blade dynamics because, in combination with aerodynamic forces, they can 
lead to flap-lag instability (refs. 1, 2). Other portions of the equations are doubly underlined; these 
are the elastic coupling terms which can be strongly stabilizing for the flap-lag motion (ref. 2). 


Aerodynamic Forces 

The aerodynamic forces per unit length denoted by Ly and L z are the components perpendicular 
and parallel to the plane of rotation, and are related to the lift and drag forces of the rotor airfoil by 
the following. 


Ly = -L sin <p - D cos <p 
L z = L cos p - D sin p 


(34) 

(35) 


The angle defines the angle of the resultant instantaneous velocity vector V with respect to 
the plane of rotation. These parameters are illustrated in figure 3. The angle p is defined by the two 
orthogonal components of V 


U P 

p = tan 1 — 

u T 

(36) 

The equations for lift and drag based on quasi- 
steady, strip theory aerodynamics are 

L = - p acV 2 a. 
2 °° 

(37) 

D = - p c V 2 c j 
2 °° d o 

(38) 




Figure 3.- Orientation of fluid velocity 
components and aerodynamic forces. 
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where the local section angle of attack is equal to 

a=8 -ip 


( 39 ) 


Note that the local blade pitch angle 6 is identical to the pitch angle imposed at the blade root 
since torsional elastic deformations are not treated. However, small geometric rotations of the blade 
section chord line due to combined flap and lead-lag bending are neglected. (Although these rotations 
are very small, they can become significant when blade pre-cone is present.) Combining equations 
(34) through (39) and using small angle approximations yields the following equations for the 
aerodynamic forces: 



(40) 


(41) 


(42) 


dv 

U-r = Six + — 
T dt 


(43) 


where v i is the induced velocity of the rotor. A reasonably accurate equation for the induced velocity 
of a hovering rotor can be obtained from blade-element momentum theory (ref. 4). Approximating 
the induced velocity as constant over the entire area of the disc, equal to the value at £ = 3/4, we 
obtain: 


SIR 



(44) 


After noting that ( UpjUj ) 2 « 1 and c^Ja « 1, the dimensionless relations for the aerody- 
namic forces become 




¥ + 


( W < +2{C -?) 


dv 

9 T 


+ (6Z- 2v t ) 




d? -^ + ( 20 $-^) 


dv 

dr 



(45) 

(46) 
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The equilibrium and perturbation components of the aerodynamic forces to be used in equa- 
tions (28) and (33) can be obtained from equations (45) and (46) and are 


■>° 6 




(47a) 


L zo = 


(47b) 


Al y-l 


(*i**t) 


.... . 3Av _ 3 Aw 

- ) 17 


3 r 


(48a) 


7 f" 3Av 3Aw*l 

— + £— J 


(48b) 


Final Equations 

Combining equations (28) and (47) yields the final equilibrium equations 


dv 0 

“d? 


/l - ^ 2 \d 2 v 0 

\ 2 ) d£ 2 


+ a 2 

2 A 2 d| 4 


d 4 v 0 , -r d 4 w 0 _ _ 7 1- 


+ Aj 


2 - "0 


d£ 


- V n = 




(49a) 


S 


dw r 


(r) 


d 2 w 0 d 4 v 0 2 d 4 w 0 _ 7 

d £ 4 1 d £ 4 6 v “ s s "' v 


t2 +A 12 


m 2 -tVi)-P pc i; 


dS \ 2 / d£ 

The final perturbation equations obtained from equations (33) and (48) are 
3AE 




/l-S 2 \a 2 Av 2df 0 3Av 2d 2 v 0 f ' 3A7 _ 2 3 4 Av 3 4 Aw 

(2 j ~W + “d i Yt w l 3r ^ + 2 3S 4 + 1 2 3S 4 


(49b) 


4 Aw 3 2 Av 
+ 


3£ 4 3 r 2 


A _ „„ 3Aw 3 /•Wdi'o 3Av dw 0 3Aw\ 7 

pc 3r a t J \ d£ 1 ' 


3£ d£ 3? 


/ _ ^dn\ 3Av 

h +2 ^h7 


_ 3Aw . 

+ (0S-2v / ) — |= 0 


3r 


(50a) 
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3 Aw (\ - 3 2 Aw ^ 2dw 0 3Av 2d 2 w 0 /* 3Av , 

* UT “ y 2 " j 3| + “dT ” d£ 2 J 17 


3Av _ 3 4 Av _ 3 4 Aw 3 2 Aw 
— d? + A 12 -gpr + Ai -gp- + -^r- 


3Av 7 


3Av 3Aw 


+ 2 ^ ,7+j o^e =0 (50b) 


Solution of Equations of Motion 

The flap-lag equations that have been developed can be further simplified using Galerkin’s 
method. The nonhomogeneous equilibrium equations are linear, variable coefficient, ordinary 
differential equations for v 0 (£) and w 0 (£). The homogeneous perturbation equations are linear, 
variable coefficient integro-partial differential equations for Av(£,r) and Aw(£,r). The coefficients 
of the perturbation equations are also functions of the equilibrium displacements v 0 , w 0 as a con- 
sequence of the nonlinearity of the original equations of motion. For Galerkin’s method we assume 
an approximate series solution in terms of the elastic bending modes i pj{%) 


j. v 

-Z 


- Z V' <{) 


ai7 = Z A ¥ T W (> 


i v 

Aw = ]Ta 


These series expressions are substituted in equations (49) and (50), which are multiplied in turn 
by each bending mode function t^(£) and integrated along the length of the beam from £ = 0 to 1. 
This procedure reduces the equilibrium equations to 2N linear algebraic equations for ? 0 . and 
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/L , and the perturbation equations to 2 N linear constant coefficient ordinary differential equations 

/ 

for A f -(|) and A/ 3 y (£). The integrals involving various combinations and derivatives of <^(f) are the 
Galerkin integrals. The equation for the bending modes used to evaluate these integrals is 


<p ; -(£) = cosh(c op - cos(a> ; £) - /ytsinhfojyl) - sin(cu^)] (52) 

This is the equation for the modes of a nonrotating uniform cantilever beam, and the constants 
co - and fj available in reference 5 are given in table 1. The Galerkin integrals are evaluated exactly, 
wnere possible, according to the tabulated integrals of reference 6. Otherwise they are evaluated by 
numerical integration. 

TABLE 1 .- CONSTANTS FOR NONROTATING CANTILEVERED • 


BEAM BENDING MODES 


j 

CO- 

i 

f J 

1 

1.875104068712 

0.7340955137589 

2 

4.694091132974 

1.018467318759 

3 

7.854757438238 

.9992244965174 

4 

10.99554073488 

1.000033553252 

5 

14.13716839105 

.9999985501087 


The equilibrium and perturbation equations are then 



J= l 


(53a) 


N 

/=' 





(53b) 


N 

E 

7=1 


Af y-(5fy) + Afy 


*V|( 2 V V***) 


+ A0 # 


~ ^Ppc^ij ~ 27 ij ~z (QEij 2v Ay) 


+ Afyl/)^ + (A 2 2 co ; - 4 - 1 )5 /y ] + A/3 y -(Aj 2 co y 4 5 /y -) 


= 0 ( 54a ) 
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N 

£ 


==1 


r 


A W + A *> 


2fipc^ + 2Tjt+l (yfig-USq) 


+ A/3,- 



I/) + A ^/ A l 2 w y Sy. 


+ A^.(Z) // + A 1 2 co j 4 5 // ) 


= 0 1 = 1 , 2 ,. 


,iV 


where the Galerkin integrals are 



D U 





V?Pj" d£ = 


o 



'p/'Pf d| 


F ijk - J“4 J 'PfPf'Pk" d £ + / ^/V/c dt l - J 'Pi'Pj'Pk'" d£ 

* o o o 

(note J" ^ J */* k 'V = -± f <# = F kji ) 

0 0 * 0 


(54b) 


(55a) 


(55b) 


(55c) 


(55d) 


(55e) 


(550 
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and where 



(Fikj Fjki)$ok 


(56a) 



'(56b) 


These equations are solved using standard methods. The equilibrium displacements are used to 
determine the coefficients of the perturbation equations, which are then solved for eigenvalues and 
thus the stability characteristics of interest. 


RESULTS AND DISCUSSION 


Uncoupled Natural Frequencies 

When the aerodynamic terms in the rotor blade equations are eliminated, the resulting perturba- 
tion solutions provide the rotating natural frequencies of the rotor blade motion. For untwisted 
blades at zero pitch angle, the equations for the flap and lead-lag degrees of freedom are independent, 
and the solutions yield the uncoupled natural frequencies. 

The first two flap and first two lead-lag uncoupled rotating frequencies of typical hingeless 
rotor blades are given in figure 4. Two examples are shown, one a soft inplane configuration and the 
other stiff inplane. The terms soft or stiff inplane refer, respectively, to rotor blades with first rotat- 
ing lead-lag frequencies less than or greater than the nominal rotational frequency of the rotor. 
These results, using three bending modes for both the flap and lead-lag degrees of freedom, illustrate 
typical variations of the rotating frequencies with the rotor angular velocity £2. The frequencies and 
angular velocity are made dimensionless by a nominal angular velocity 12 0 , which can be interpreted 
as the normal rotational velocity of the rotor. The dimensionless bending stiffness of the rotor 
blade is inversely proportional to the rotational speed, K = K 0 /(n/£l 0 ) 2 where K Q = EA/m£l 0 2 R 2 . 

The convergence of the solution with the number of bending modes is shown in figure 5. Here, 
the rotating frequencies are shown as a function of the nonrotating frequencies. They are both made 
dimensionless by the rotor angular velocity. The use of this nondimensional parameter is a means of 
distinguishing between different rotor blade configurations. For a given rotor speed, then, the non- 
rotating frequency is a measure of the bending stiffness of the rotor blade. The first flap and lead-lag. 
frequencies are of primary interest here and the results show that convergence is quite good for rotor 
blades of high bending stiffness (or low angular velocity) but becomes poor for low stiffness (or high 
angular velocity). To give proper perspective to these results, we must consider them in the light of 
typical blade frequency values. For the lead-lag degree of freedom, a typical stiff inplane hingeless 
rotor blade would have 65^ ~ 1.5, and a single bending mode would accurately define the rotating 
lead-lag frequency. For a soft inplane configuration, 0.7, two bending modes would be required. 
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Figure 4 — Variation of the first two uncoupled rotating natural frequencies with rotor angular velocity 6 = 


Dimensionless rotating natural flap frequency, 





For the flapping degree of freedom, a typical first bending mode frequency would be co^ ~ 1.1, 
which would require at least two bending modes. This is particularly important for hingeless rotors 
because the flapping response characteristics can be sensitive to small changes in the first flap mode 
frequency. 

It should be noted that the convergence characteristics in figure 5 are directly related to the use 
of the simple nonrotating uniform cantilever beam mode shapes and can be viewed as a measure of 
the usefulness of approximate modes. However, for use in the stability analyses reported herein, the 
use of nonrotating beam modes simplifies the parametric studies because different mode shapes are 
not required for different rotor blade stiffness values. 

Also shown for comparison in figure 5 are the rotating frequencies obtained using the approxi- 
mate centrally hinged, spring restrained, rigid blade representation. The uncoupled flap and lead-lag 
natural frequencies are given by the following relations (ref. 2): 


“(j = vH* <57a) 

“{■“"Saw (57b) 


In effect, these frequencies are a special result for the case of linear mode shapes for the flap and 
lead-lag bending deflections. These modes might be expected to yield the correct first rotating fre- 
quencies in the limiting case of vanishing blade bending stiffness because the mode shape approaches 
a straight line. It must be noted, however, that because of the condition of zero slope at the root of 
the cantilevered blade a linear mode shape is invalid at that point. 


Stability Characteristics 

The dynamic stability characteristics of the flap and lead-lag bending oscillations of elastic 
hingeless rotor blades without pre-cone in a hovering flight condition is shown in figure 6. The results 
are presented in locus of root form for both the first flap and lead-lag modes as the blade pitch angle 
is increased from zero. Loci for several configurations having various lead-lag stiffnesses are presented. 
As discussed above, only the soft and stiff inplane configurations are representative of actual hinge- 
less rotors. The other parameters are chosen to be representative. 

The flap mode is typically well damped. The lead-lag damping is substantially increased for 
higher pitch angles with no indication of instability. This is consistent with earlier findings of refer- 
ence 2 for the rigid blade model with the appropriate elastic coupling characteristics. In effect, the 
elastic flap-lag coupling introduces the large flap mode damping into the lead-lag mode, thereby 
augmenting the small inherent lead-lag damping and nullifying the destabilizing inertial and aerody- 
namic flap-lag coupling effects. 
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Figure 6.- Locus of roots with increasing pitch angle, y = 5.0, 
ct = 0.05, 0.6, c^ = 0.01, N= 3. 


A comparison between the results of the elastic blade and the rigid blade model will now be 
made. A schematic picture of the approximate centrally hinged, spring restrained, rigid blade is given 
in figure 7. The orientation of the flap and lead-lag restraint springs (parallel and perpendicular to 
the blade chord line, respectively) simulates the elastic coupling characteristics of the actual elastic 
blade. The spring stiffnesses are chosen so that the uncoupled rotating flap and lead-lag natural 
frequencies match the corresponding first mode rotating natural frequencies of the elastic blade. 
The roots for the rigid blade configuration are shown in figure 8, and they exhibit behavior nearly 
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Figure 7.— Schematic of approximate centrally hinged, spring restrained, rigid blade 
representation with springs oriented to match the elastic coupling of elastic 
rotor blades. 

identical to the elastic blade results discussed above. The differences are a reduction in flap mode 
damping and a small reduction in the stiff inplane lead-lag damping for the rigid blade representa- 
tion. The significance of the close agreement is that the mode shape does not play a large role in 
the flap-lag stability characteristics, at least in the hover condition, and the simple rigid blade model 
should be quite adequate for approximate stability analyses if the equivalent spring constants are 
chosen as discussed above. 

Although the results just presented confirm the relatively small effect of mode shape on stability, 
it is of interest to examine the sensitivity of the elastic blade stability to the number of flap and lead- 
lag bending modes included in the solution. The roots in figure 6 were computed using three bending 
modes for each degree of freedom; we will now investigate the effect of using only one or two modes. 
Since the first lead-lag mode damping is of most interest, we will consider only the real part of this 
root. Betore presenting the results it should be noted that there are two ways of making the compar- 
isons. This depends on whether the rotating or nonrotating first mode bending frequency is consid- 
ered to be the independent parameter defining the rotor blade configuration. If one uses the rotating 
frequency as the independent parameter, as for the rigid blade comparison, the nonrotating fre- 
quency, or blade stiffness parameter K, will depend on the number of bending modes used. The 
particular variation can be obtained from figure 4. In this case the effect of the number of bending 
modes on lead-lag damping will be primarily due to coupling between the various bending modes. 
Alternatively, if the blade stiffness parameter K or the nonrotating frequency is chosen as the inde- 
pendent parameter, variations in the number of bending modes will directly influence the rotating 
first mode frequency. This would introduce an additional effect, since the lead-lag damping is influ- 
enced by the proximity of the flap and lead-lag rotating frequencies. 

The comparisons of first-lead lag mode damping are given in table 2 for a soft and a stiff inplane 
rotor blade configuration operating at a pitch angle of 0.3 rad. The sensitivity of damping to the 
number of modes is quite small but certain trends are discernible. As expected, holding the rotating 
natural frequencies constant minimized the effect of modal interaction, in contrast to holding the 
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TABLE 2.- EFFECT OF NUMBER OF FLAP AND LEAD-LAG BENDING MODES 
ON THE FIRST LEAD-LAG MODE DAMPING AT 9= 0.3, Op = 0.4. 


Soft inplane configuration, 

aJs- = 0.6 

>NR 

Stiff inplane configuration, 

c3 > =1.4 

f NR 

N 

-a^/n 

N 


3 

0.01689 

3 

0.02342 


ojf- const 
*NR 

a>£ const 


c Of- const 

^NR 

ajg const 

2 

0.01663 

0.01694 

2 

0.02352 

0.02344 

1 

.01567 

.01714 

1 

.02390 

.02363 

Rigid 


.01713 

Rigid 

... 

.02253 



Figure 8.— Locus of roots with increasing pitch angle for rigid 
blade representation. Rotating natural frequencies are cho- 
sen to match the corresponding results in figure 6. 


nonrotating frequencies constant. Further- 
more, the effect of number of modes was 
more evident for the soft inplane configu- 
ration, which is consistent with the earlier 
results. This can be attributed to the increas- 
ing disparity between nonrotating and rotat- 
ing mode shapes as the stiffness is reduced. 
For completeness, the table includes results 
presented graphically in figure 8 for the 
approximate rigid blade. It is interesting to 
note that the rigid blade is as accurate as 
the single mode solution for the soft inplane 
configuration, but underestimates the damp- 
ing for the stiff configuration. Results for 
the rigid blade for the case of constant non- 
rotating frequency are not presented because 
of the unrealistic frequencies of this approxi- 
mate representation. 

These results indicate that the basic 
character of flap-lag stability is little changed 
by using the approximate rigid blade (with 
the proper rotating natural frequencies) or 
a single mode representation for the elastic 
bending deflections. A small degradation in 
the accuracy of the lead-lag damping can 
result from the rigid blade approximation, 
but a single bending mode gives quite 
reasonable accuracy. 

The nonlinear Coriolis and centrifugal 
terms in the equations discussed above were 
considered essential to properly describe 
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the flap-lag stability characteristics of elastic blades. Although the above results indicate that the 
destabilizing aerodynamic and inertial coupling is more than compensated for by the elastic coupling 
of a rotating uniform beam, the destabilizing nonlinear terms are also of interest for their influence 
on the lead-lag damping. This is shown in figure 9, again for two rotor blade configurations. The 
nonlinear terms do act to destabilize the lead-lag degree of freedom and their effect is substantial. 
For rotating beam configurations with radially nonuniform stiffness distributions, reference 2 implies 
that the elastic coupling may not, in general, prevent flap-lag instability, and the nonlinear inertial 
terms will then determine the rotor blade stability characteristics. This problem, however, is beyond 
the scope of the present report. 



(a) Soft inplane configuration, cjj. ^ = 0.6. 


Figure 9.— Effect of the nonlinear inertial terms on the lead-lag mode damping, 
7^5.0, a = 0.Q5,^ NR : 


= 0.6,c d = 0.01, N = 3. 
u o 
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Rotor blade pre-cone is an important design parameter for hingeless rotors and previous analyses 
(refs. 1 and 2) with the approximate rigid blade representation indicate that it has a substantial effect 
on flap-lag stability. Some representative results for the lead-lag damping of the elastic blade are 
shown in figures 10(a) and 10(b) for the soft and stiff inplane configurations, respectively. In general, 
either positive or negative pre-cone is stabilizing, except for small positive increments with the stiff 
inplane configuration. 



(a) Soft inplane configuration, ^ =0.6. 

Figure 10.— Effect of pre-cone on the lead-lag mode damping, 7 = 5.0, 
a =0.05, (5^ = 0.6, c d =0.01, N= 1 , 8 = 0.1. 
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(b) Stiff inplane configuration, =1.5. 

Figure 10.— Concluded. 

The effects of rotary inertia terms in the rotating beam equations were determined to be of 
nearly negligible magnitude. This was checked with the numerical results by computing the lead-lag 
mode damping both with and without the rotary inertia terms. The results are given in table 3 and 
confirm the negligible effect on flap-lag stability. 


TABLE 3.- EFFECT OF ROTARY INERTIA ON THE FIRST LEAD-LAG MODE DAMPING, = 0.6, 

%R =1-5,7= 5.0, c do = 0.01 , a = 0.05 , 6 = 0.3, N = 1 



No rotary inertia 

With rotary inertia 

0.021199 

0.021190 
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CONCLUDING REMARKS 


1 . Conventional linear equations of motion cannot be used to determine the stability of coupled 
flap and lead-lag degrees of freedom of cantilevered, torsionally rigid, elastic rotor blades. An addi- 
tional equation for the radial degree of freedom and a nonlinear strain-displacement relation are 
required. The basic nonlinear flap-lag equations can be linearized about a suitable equilibrium condi- 
tion and still retain the effect of the important inertial flap-lag coupling terms. The calculation of 
lead-lag damping was shown to be very much a function of these terms. 

2. Flap-lag oscillations in hover for untwisted blades of uniform mass and stiffness were found 
to be stable over a wide range of parameters. It was determined that the effects of elastic flap-lag 
coupling introduced sufficient flap mode aerodynamic damping into the weakly damped lead-lag 
mode to overcome the destabilizing inertial coupling. Previous results (ref. 2) for an approximate, 
centrally hinged, spring restrained, rigid blade representation, however, indicate that this effect can 
be reduced for elastic blades with twist and/or nonuniform mass and stiffness. 

3. The first mode frequency and damping for the elastic blade motions were found to be nearly 
identical to results obtained with the rigid blade representation provided that the rotating natural 
frequencies are equal. This finding enhances the validity of the simpler rigid blade model for investi- 
gating the flap-lag stability characteristics of hingeless rotor blades in hover. 

4. The effect of the number of elastic bending modes on the stability is generally small, espe- 
cially where the rotor blade stiffness parameter is defined by the rotating natural frequency. Very 
accurate results for the lead-lag damping can be obtained with two modes for both the flap and 
lead-lag degrees of freedom. 

5. For practical hingeless rotor blade configurations, the fundamental uncoupled flap and lead- 
lag rotating natural frequencies can be determined with acceptable accuracy by retaining the first 
two nonrotating bending modes for each degree of freedom. 

6. Rotor blade pre-cone is generally stabilizing, although in some cases it is moderately 
destabilizing. 

7. Rotary inertia effects are negligible for practical rotor blade configurations. 
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